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Abstract 



The Bethe equations, arising in description of the spectrum of the dilatation operator 
for the su(2) sector of the JV = 4 supersymmetric Yang-Mills theory, are considered in the 
anti -ferromagnetic regime. These equations are deformation of those for the Heisenberg 
XXX magnet. It is proven that in the thermodynamic limit roots of the deformed 
equations group into strings. It is proven that the corresponding Yang's action is convex, 
which implies uniqueness of solution for centers of the strings. The state formed of strings 
of length (2n+l) is considered and the density of their distribution is found. It is shown 
that the energy of such a state decreases as n grows. It is observed that non-analyticity 
of the left hand side of the Bethe equations leads to an additional contribution to the 
density and energy of strings of even length. Whence it is concluded that the structure of 
the anti-ferromagnetic vacuum is determined by the behaviour of exponential corrections 
to string solutions in the thermodynamic limit and possibly involves strings of length 2. 

1 Introduction 

Integrable models, in particular, spin chains, appear in several problems of high energy 
physics as effective models of interaction. For example, the Hamiltonian of the XXX spin 
chain with the non-compact representation of spin s= — 1 arises in description of scattering 
of hadrons at high energies [H [2] , and also in description of mixing of composite operators 
under renormalization in QCD [3]. 

Mixing of composite operators under renormalization in the super-symmetric Yang- 
Mills theory also gives rise to an XXX-chain but with a compact representation. In this 
theory, one considers locally invariant operators of the form 



where Z and W are two complex scalar fields from the supermultiplet. The conformal 
dimensions A of this operators comprise the spectrum of the dilatation operator D. It 
is convenient to describe mixing of operators ([1]) under renormalization with the help of 
an analogy with the quantum spin chain of length L = J\ + J2, where each occurrence 
of Z is represented by a spin up, and each occurrence of W is is represented by a spin 
down. For example, the state ZZZWWZWZ corresponds to the following spin chain 
TTT11T1T I)- An important observation made in [3] was that, in the su(2)-sector of the 
theory in the one-loop approximation (i.e., in the first order in A = 5y M Y, where g YM is 
the Yang-Mills coupling constant, and N is the number of colours) the dilatation operator 
D can be expressed via the XXX Hamiltonian of spin s=^ for the described above chain: 



O = tr{Z Jl W j2 + permutations), 





(2) 



Therefore, determining the spectrum of D in this approximation is reduced to investigating 
the Bethe equations for the Heisenberg magnet (see El Ej ) , 

Uj +j/2 \ L _ yr Uj - u k +i 

Uj-i/2) l^Uj-Uk-i 1 ' ,j 

where Uj are the rapidities of elementary excitations. 

The high-loop corrections to ([2|) were found in articles that followed [4 J , and there it was 
shown that the corresponding expressions are also integrable Hamiltonians for the spin 
chain (those that include interaction between several nearest sites). In these approxima- 
tions, the spectrum of D is determined not by equations ([3]) but by their "deformations" 
which explicitely contain the parameter A in the left-hand side. Assuming that integra- 
bility of D takes place in all orders, Beisert, Dippel, and Staudacher [8] argued that the 
exact Bethe equations determining the spectrum of D look like following 

x(uj + i/2) \ L _TT u j ~ u k + i 

x(uj - i/2) J Uj -u k -i' 

where 

x(u) = - (u + y/u 2 - x 2 ) X 2 = ^ ' ( 5 ) 

For the spin chain of length L, the spectrum of the Hamiltonian is bounded by the 
energies of the ferromagnetic and the antiferromagnetic vacua. Therefore, the spectrum 
of the operator D, i.e. the dimensions of operators ([!]), belongs to the interval 

L < A sC A max . (6) 

The lower bound here follows obviously from the structure of the antiferromagnetic vac- 
uum (all spins up), whereas determining the upper bound requires a non-trivial evaluation 
of the energy of the antiferromagnetic vacuum. The value of A max in the thermodynamic 
limit was found in [9] and [10] by means of the standard technique of passing to the limit 
L — > oo in the Bethe equations. In this procedure it was assumed that, like for equations 
([3]), the antiferromagnetic vacuum is formed of strings of length 1, i.e., real roots of equa- 
tions However, the example of the XXX spin chain of higher spin [11,12] shows that 
the antiferromagnetic vacuum for Bethe equations whose left-hand side differs from (J3j) 
can have other structure, for instance, it can be filled by strings of greater length. 

The aim of the present work is to prove existence of solutions corresponding to strings 
of length greater than 1 for the equations (JIJ) in the thermodynamic limit and to check 
validity of the assumption about the structure of the antiferromagnetic vacuum. 



2 Existence of string solutions 

To make from x(u) an analytic single-valued function, we fix in C the cut [— XiX] an d 
define x(u) on C/[— x, x] as follows: 

x(u) = i(Vnei ei + ^e^) 2 , (7) 

where r\, T2 £ M+ and 9\, 62 6 ]— vr, ir] are determined from the relations u = x + fi& tdl = 
— x+f2e lS2 . Notice that the signs of the imaginary parts of u and a; (it) coincide, that is x(u) 
maps a point from the upper/lower half-plane to the upper/lower half-plane, respectively. 
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Let us prove that the following relations 



> 1, for Imu > 0, 

= 1, for lmu = 0, (8) 
< 1, for Imu < 0. 

hold for every 5 > 0. 

Let us denote s = Re(u), t = Im(-u), a = Re(x(ii)), b = Im(x(-u)). It follows from (7) 
that |x(n)| is a function continuous in u (in particular, even when u crosses the cut), and 
\x(u)\ = \x(u)\. Therefore, with s being fixed, the function \x(s + it)\ is continuous and 
symmetric in t. We will prove that this function is convex, of which (8) will then be an 
obvious consequence. 

2 

It follows from (5) that the function inverse to x(u) is given by u(x) = x + Whence 
it is easy to derive the following relations 

Let dt denote the partial derivative w.r.t t (i.e. dts = 0). Applying it to (9) and solving 
the system of equations for d t a and d t b, we find 

8 t a= 2 -^, d t b=± ((a 2 + b 2 ) 2 + X 2 (b 2 -a 2 )), (10) 
where D = ((% - a) 2 + b 2 ){(\ + a) 2 + b 2 ). Whence we obtain 

d t \x(u)\ 2 = 2^±^(a 2 + b 2 + X 2 )b. (11) 

Due to the remark made after equation (7), t and b are of the same sign. Therefore, 
expression (11) is positive/negative in the upper/lower half-plane, respectively. Hence 
|x(m)| is convex in Im(n), which completes the proof of (8). 

Relations (8) allow us to adapt the analysis of complex roots of equations (3) (see [5,7]) 
to the case of equations (4). Namely, it follows from (8) that, in the L — > oo limit, the 
absolute value of the l.h.s. of (4) tends to oo when Im(uj) > and to when Im(uj) < 0. 
This implies that r.h.s. has, respectively, a pole or a zero, i.e. there must exist also the 
root Uj — i in the first case and the root Uj + i in the second case. Thus, like in case of 
the XXX magnet, roots of equations (4) group in the thermodynamic limit into "strings" 
which are complexes of the form Uj >m = Uj + im, where Uj 6 M and 2m G Z. 



x(u + iS) 
x(u — iS) 



3 Strings of odd length 

3.1 Equation for the centers of strings 

Now we will investigate what state has the maximal energy in the case when the vacuum 
is filled with strings of length 2n+ 1, where n is integer. The number of these strings v n is 
fixed by the condition (2n + l)v n = L/2. Following [5-7], we multiply Bethe equations (4) 
along a string of length 2n + 1. Since the right-hand side of these equations is the same as 
in the "undeformed" Bethe equations, strings will have the same form, i.e. Uj = n™ + im, 
m £ Z. Further, considering the thermodynamic limit (L — > oo), we obtain that the 
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centers of strings are arranged along the real axis with some density which satisfies certain 
integral equation. Having found this density, one can compute the energy of the ground 
state (see [5,6]). 

Thus, we obtain the following equation for the centers of strings: 

1 z(< + (2n + l)|, 
L i n A = + V ® n , n {u] - u n k ), (12) 

2 x (yn-{2n + l)^ ti 



where 



2n-l 



* n n (u) = arctan 1-2 > arctan . (13) 

' V ' 2n + l ^ m + 1 y J 

3.2 Yang's action 

To prove the uniqueness of a solution to equations (12) for a given set of integer numbers 
Q™, we will use the Yang's action. As in the case of the XXX magnet [7], there exists a 
functional S (called Yang's action) such that equations (12) are the conditions of its ex- 
tremum, d Ua S = 0. Let us consider the quadratic form for the matrix of second derivatives 
of S: 

d 2 S * r V-~ , x(u a + i(2n + l)) ^ 1 

|> fltfC' = * L ? ^ - i (*» + D) + 5 K - + i Vf>) ' 

(14) 

where v a £ M. It is obvious that the second term is always positive. Let us prove the 
positivity of the first one: 

1 x(s + it) a bd s a — ad s b adta + bdtb dAxis + it)\ 2 
-d s \og— -r = d s arctan - = = = . (15) 

2 x(s — it) b (a z + o" 1 ) (cr + o 2 J 2(o z + b z ) 

Here we used the same notation as in Section 2 and applied the Cauchy equations for 
derivatives of an analytic function. Relation (11) shows that (15) is positive for b > 0. 
Therefore the quadratic form (14) is positive definite and, consequently, the action S has 
a unique minimum. 

3.3 Thermodynamic limit 

Let us go to the thermodynamic limit now. Taking L — > oo and differentiating (12) with 
respect to u, we obtain for the left-hand side: 

l.h.s. = -( ; 1 1 V (16) 

Now we introduce the root density p(u): 

= > (17) 



dq) 



q=q(u) 
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where q(u) = Qj/L. p(u) plays the role of density of numbers q(u) on the interval du. 
Having introduced this density, we can rewrite the l.h.s. of (12) as follows 



2n-l 

In + 1 „ v-^ m + 1 



r.h.s. = TTp(u) / dpp(ri) - — -tt + 2 V 



n (m+l) 2 + (n-M) 2 

m=0 



(18) 



This integral equation can be solved by means of the Fourier transform. To this end we 
compute first 



f°° du ^ = e(-kl) e"l fc/ l / ^ exp (ik (q + ^X) 



y/{u + il) 2 - X 2 

,-2-k 

= sign (-l)9(-kl)e-W / id^e ikxcoSip = 2ni sign (-1)9 (-kl)e~^ J ( X k), (19) 

Jo 



where Jo(k) is the Bessel function of the first kind. Thus, after the Fourier transform, the 
left-hand side of our integral equation acquires the following form 

F[l.h.s] = Tre-^ kl Mxk) ■ (20) 
On the r.h.s. the Fourier integral is divided into the following terms 

oo oo oo 

J dpp(p) I * iku A2 + ^_ fi)2 du = ir J dpp(p)e^ k e^ A = 7Te'\ k \ A p(k), (21) 

— oo — oo — oo 

where p(k) stands for the Fourier transform of the density p(u). These terms can be 
summed up: 

/ 2™- 1 \ /o -2n|fe| „-(2n+l)|fc|\ 

np(k) I 2 e-l fc K m+1 ) +e"l fc K 2 " +1 )J = irp(k) j . (22 ) 

As a result, we obtain the following expression for the Fourier transformed density p(k): 



c 



fc l-l J (xfc)tanh^ 



p(k) - Jo(xfc) e _^ |fc|(i + e | fc|)(e(aB+1)W _ 1} " 2ginh (n + i) ]k] ■ ( 23 ) 
This yields the solution of the integral equation, 

OO |i| 

, x If Jo(xfe)tanhLI 

p(u) = — / t i"t"~~ ~ e dk. (24) 

PV ' 2vr J 2sinh(n + ±) \k\ V ' 

— oo 

For n = (i.e. strings of length 1) this expression coincides with the expression for the 
density that was obtained in [9, 10]. 

3.4 Energy of strings of odd length 

Now we will compute the energy (or, equivalently, the maximal dimension of the dilatation 
operator A max ) for the state filled with strings of length 2n + 1, where n is integer. The 
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dispersion in the considered theory differs from that of the XXX magnet, the corresponding 
energy density is given by [9] 



1 + 



00 

^2 / du P( u ) 



1 



1 



(„ + 2"+Ii) x ( 



U 



2n±l, 



(25) 



Substituting here the expression (24) for the density, we obtain: 



max _ 1 I" dk J (xfc)Ji(xfc)tanh| 

L 7T / k e ( 2n+1 ) fc -l 



it J k 




(26) 



Here we used the following relation: 

X*_ I i i 

2 



\xd x 



X' 



2i \x(u + il) x(u — il) 



(27) 



Below we provide a plot which shows how the second term in (26) depends on n for a fixed 
value of A. 




It is apparent that the maximum of the integral (and thus the maximum of A max ) is at- 
tained at n = 0. In the Appendix we give a strict proof that A max decreases monotonously 
as n grows. 

Thus, we conclude that the state with maximal energy in the sector of strings of odd 
length indeed corresponds to strings of length 1, as it was assumed in [9, 10]. However, in 
the next section we will show that the real anti-ferromagnetic vacuum might have a more 
complicated structure and possibly corresponds to strings of length 2. 



4 Strings of even length 

Unlike the absolute value |x(u)|, the phase of x(u) is not a continuous function. Its value 
changes by a finite amount when u crosses the cut. Using (7) it is easy to derive that for 
u G [— Xi x] we have the following relation: 

, x(u — ie) \J y 2 — u 2 , . 

lim In —. '- = 2i v e arctan -^A ; (28) 

e^O X\U + 16) U 

where v e = — 1 if e tends to the zero from the right and v e = 1 if e tends to the zero from 
the left. 
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For a chain of a large but finite length L, the roots of Bethe equations group into 
strings only up to exponential deviations: Uj >m = Uj + i(m + e m ), where m is integer 
or half-integer and e m = 0(e~ WmL ), uj m > 0. The total energy of a string is given by 

£j = | Yl( j^. 1 ~~ x ^ u . 1 • Computing this expression for a string of even length 
and taking relation (28) into account, we observe that the roots 

uj-i^+e) (29) 

give additional contributions, 

_ if 1 1 1 1 \ 

fj ' " So 2 U(«j-*<0 ~ x( Uj -+ie) + + X (u 3 +^M) ~ x ( Uj -2™±±ij)- (30) 

It is important to remark here that, for the regime corresponding to v e = —1, the first 
two terms in (30) give a negative contribution to £j. Furthermore, it can be shown that 
in this regime £j is positive not everywhere on the real axis. As a consequence, the anti- 
ferromagnet vacuum in this case has a more complicated structure — it has to be filled 
with strings only in the intervals where £ j > 0. In the present work we will consider only 
the regime corresponding to v e = 1. Let us remark that in this case the proof of convexity 
of the Yang's action given in Section 3.2 remains valid. 

Taking product of the left hand sides of the Bethe equations (4) along a string of even 
length and taking into account the additional contributions due to the roots (29), we find 
that the l.h.s. of the equation for the centers of strings looks like following 

h M, rii )W,(», + |( ! » + i)A 1 

e/Q \x{ Uj + ie)) y x { Uj - i(2n + 1)) J 

Taking logarithm of this expression, we obtain equations (12) but with an additional term 
on the left hand side. Differentiating (28) w.r.t. u, then making the Fourier transformation, 



d u lim In X ^ 



/x e iku ^ u 
, = = = -2ttz v e J ( X k), (32) 



e^o x{u + ie)\ ' i_ x _ 

and comparing with formula (19), we see that, for strings of even length in the regime 
v t = 1, equation (20) acquires the following form 

F[l.h.s] = 7r(e-^l fe l + l)Mxk). (33) 

Repeating now the computation of the Fourier image of the density as was done in Sec- 
tion 3.3, we find 

— oo 

Formula (33) shows that the additional contributions from the roots (29) yield the same 
effect as if the vacuum were filled with mixture of strings of length (2n+l) and strings 
of length (which formally corresponds to n = -5). In this context we remark that 
expressions of the type (31) appear in the left hand sides of Bethe equations for chains 
with alternating spins (see, e.g. [13, 14]). 
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Formula (30) implies that, in order to compute the total energy of the chain, we have 
to replace (25) with 

OO \ 

A _ ix 2 . f r ,( 1 1 1 1 \ 

— OO 

(35) 

Substituting here (34) and making the same computation in the Section 3.4, we obtain 

OO 

A f dk T , ,, T . ,, , k , (2n+l)fc , , 

- = 1 + 2 X / — Mxk)Ji(xk) tanh - coth V - ; . (36) 

o 

Using the method described in the Appendix, one can show that for n > \ expression (36) 
decreases monotonously as n grows. The maximal value of (36) at n = \ is given by 

OO 

^p = l + 2 X f fj (xk)Ji(xk). (37) 
o 

In order to compare this expression with the value of *g ax corresponding to strings of 
length 1, one can subtract the value of (26) for n = from (37). The resulting expression 
has the form of the integral in (38) with monotonously decreasing in k function f(k). As 
shown in the Appendix, such an integral is positive. 

Thus, filling the vacuum with strings of length 2 in the regime v e = 1 yields larger 
value for Al j? ax than filling it with strings of length 1. However, what regime is indeed 
realized for string of even length in the thermodynamic limit, remains at present an open 
question. Its resolution requires quite subtle analysis of the exponential corrections e in 
(29) for L -> oo. 



Appendix 

Let us denote a = (2n + 1), x = v / A/(27r) and make in (26) a substitution k' = xk- Then 

OO 

d a (^p) = -2 J dkJ Q {k)J x {k)f{k), (38) 
o 

where 

qfc k 

1 \ / \ f e x — 1 



f(k) = — — • (39) 

Vex+1/ \ e — -1/ V e — -1/ 

In this form, it is evident that f(k) decreases monotonously as k grows for all a > 1 and 
X > 0. Using this we will show that the integral in (38) is positive. 

Let < t\ < *3 < is . . . be the ordered set of roots of Jo(t), and = to < ti < . . . be 
the ordered set of roots of J\ (t) . It follows from the relation 

d t J (t) = -J ± (t) (40) 

that t2n < ^2n+i < ^2rt+2- Taking into account that Jo(t)Ji(t) is positive on ]t2 n , *2n+i[ and 
negative on ]t2n+i, and the function f(k) is positive and monotonously decreasing 
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for all k > , we obtain the following estimate: 



00 



dk J (k) Ji(fc) f(k) = ) I dkJ (k)J 1 (k)f(k)> (41) 



n=0 J t2n 

00 '-t2n+2 < X 



"~ r ft 2n +2 -, 1 r 

> J] |/(Wl) / JlWj = g E [/(Wl) (4(*2n) " 4(W2)) 

n=0 * 2n n=0 

In the last equality we used relation (40). Now, since roots of Ji(t) are the points of local 
extrema for Jo(t), and the values of |Jo(i)| at these points form a decreasing sequence, 
we conclude that the sum on the r.h.s. of (41) and hence the initial integral are positive. 
Thus, d a ( Al ? ax ) < 0, which implies that A ? ax decreases monotonously as n grows. 
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